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Abstract 

We investigate core structures of X-ray emitting intracluster gas based on the so-called /3-model, which 
is an isothermal hydrostatic model often used in observational studies. We reconsider the j3- model and find 
that the virial temperature T vlr of a cluster may be represented better by j3Tx than Tx, where (3 is the 
parameter obtained from the X-ray surface brightness and Tx is the emission-weighted mean temperature 
of the gas. We investigate 121 clusters observed by ROSAT and ASCA and find that the luminosity- 
temperature relation Lx — (3Tx is less steep than Lx — Tx- We classify the clusters into two core-size 
groups in order to investigate their properties in detail. While in the larger core group the core radius 
is marginally proportional to the virial radius, no significant relation is found for the smaller core group. 
This may suggest that the smaller cores reflect the presence of cD galaxies, effect of radiative cooling or 
asymmetry in the surface brightness. We examine such possibilities, and find that the clusters of smaller 
cores have shorter cooling time than the Hubble time, while no significant correlation is found with cD 
or asymmetry. We carry out hydrodynamical calculations to simulate the /3-modcl, intending to see the 
behavior of the isothermal gas under the gravitational potential including the dark matter and galaxies 
with or without a central cD galaxy. Calculations show r c cx r v ; r and T v ; r ~ (3T consistently with our 
consideration to the /3-model. Also is found from calculations that the presence of a large cD galaxy may 
form a gas core ~ 40 kpc, which seems too small to account for the range of the core sizes, 40-80 kpc, of 
the smaller core group. 

Key words: galaxies: clusters: general — galaxies: evolution — X-rays: galaxies: clusters 
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1. Introduction 

Galaxy clusters are the self-gravitational systems dom- 
inated by the dark matter. The overall structure of the 
intracluster gas can be explained by the /3-model, which 
consists of a flat core part and the envelope. In the (3- 
model, the core radius and the envelope slope represented 
by (3 are important parameters to characterize the struc- 
ture of the gas. 

In simple models of structure formation by gravity, 
galaxy clusters constitute a homologous family (see e.g., 
Kaiser 1986). Clusters are self-similar in shape and their 
properties relate in predictable scaling laws. For example, 
the core radius is expected to be proportional to the virial 
radius of the cluster. 

Observations show, however, that core radii of clusters 
are not simply scaled. For 79 clusters observed by ROSAT 
and ASCA, Ota & Mitsuda (2002) find two typical core 
sizes or two core-size groups: the core radius distribution 
exhibits two distinct peaks at ~ 50 kpc and ~ 200 kpc 
for H = 70 km s" 1 Mpc" 1 (Ota & Mitsuda 2004). They 
also show that some clusters are expressed better by the 
double /3-model with two cores than the single /3-model, 
and then the radii of the larger and smaller cores fall into 
the above two core-size groups respectively. These results 
can be an important clue for understanding the formation 
and evolution of cluster cores. 

Observations also show that the mass-temperature or 



luminosity-temperature relation of clusters do not sim- 
ply scale self-similarity (e.g., Markevitch 1998; Arnaud 
& Evrard 1999; Finogucnov et al. 2001; Sanderson ct 
al. 2003). Cosmological hydrodynamics calculations sug- 
gest that various-sized clusters are formed through the 
structure formation in the universe. But such clusters are 
neither isothermal (e.g., Borgani et al. 2004) nor those 
expressed by the /3-modcl. 

As for the core structure, effects of radiative cooling can 
be a considerable factor. For instance, the cooling flows 
may increase the density of the cluster central region (see 
Fabian 1994 for a review), though recent ASCA, Chandra 
and XMM-Newton observations (e.g., Makishima et al. 
2001; Kaastraet al. 2001; Peterson et al. 2001; Lewis et al. 
2002) show several results inconsistent with the standard 
cooling flow model. Although various models have been 
proposed for thermal balance against cooling flows, (see 
e.g., Bohringer ct al. 2002; Peterson et al. 2003 references 
therein), quantitative account is not yet given successfully. 

Another factor of the double core may be the presence 
of a cD galaxy in the cluster central region (e.g., Jones 
& Forman 1984; Xu et al. 1998 and references therein). 
However, the cores associated with cDs seem small com- 
pared to the typical core size of the smaller core group; 
for example, in the Centaurus cluster, the core radius due 
to the cD is estimated to be ~ 40 kpc or less (see Ikebc 
et al. 1999). 

In the present paper, we investigate core structures of 
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the isothermal intracluster gas. The origin of the two core 
sizes is of our interest. In section 2, we examine the (3- 
modcl considering the isothermal X-ray temperature and 
the virial temperature. In section 3, we investigate 121 
clusters observed by ROSAT and ASCA and discuss the 
properties of the smaller and larger core clusters. In sec- 
tion 4, we show the results of hydrodynamical calculations 
to simulate the /3-model or the behavior of the isothermal 
gas under the cluster potential, and discuss the effect of 
a central cD galaxy. We summarize the present work in 
section 5. 

2. Isothermal /3-model 

The /3-model has been widely used in analyses of the 
X-ray emitting hot gas in clusters of galaxies. This is an 
isothermal hydrostatic model based on the approximate 
King model for galaxies in clusters (see e.g., Sarazin 1986). 
The density profile of the intracluster gas based on the f3- 
modcl is expressed as 



Pg (r)= Pg0 [l + (r/r c ) 2 }- 3P/2 , 



(1) 



where p g o is the gas density at the cluster center and r c 
the core radius. These quantities and (3 are obtained from 
the X-ray surface brightness profile for isothermal tem- 
perature Tx- The core radius r c is defined originally as a 
radius at which the projected galaxy density falls to half 
its central value (Sarazin 1986), and characterizes the gas 
distribution as well. 

The total mass enclosed within radius r is obtained as 



M{r) 



ZkTxPr (r/r c ) 2 
/j,mG 1 + (r/r c ) 2 



(2) 



where Tx is the gas temperature inferred from X-ray ob- 
servations mentioned above, fi the mean molecular weight 
and m is the hydrogen mass. Thus, the total mass in- 
tegrated up to the virial radius r v ; r of the cluster comes 
to 



M{r v 



where 



3kT x (3 1 V 



ir I Tc • 



(3) 



On the other hand, for a virialized cluster, the virial 
mass is given by 



fimG 



(4) 



where T v j r is the virial temperature. For M(r v - lr ) = M vu ., 
we have the relation between the two temperatures as 



T ■ — 



1 + x 



PT; 



(5) 



The observed values of /3x^ iT / (1 + ij ir ) = T viT /T x for dis- 
tant clusters (Ota & Mitsuda 2002; see section 3) are plot- 
ted in figure 1. 

Many of the clusters show x v i t 3> 1, and therefore 
T v i T /Tx — (3 ~ 0.6 according to equation (5). This means 
that Tx could deviate from T v ; r depending on (3, and T V j r 
is well represented by (3Tx rather than Tx when x v - 1T ^> 1. 
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Fig. 1. The ratio of the virial to isothermal temperatures is 
plotted against x v - lT . Each solid line represents the relation 
given by equation (5) for /3=0.6, 0.8, 1.0, and 1.2. The open 
circles represent /3x^ iT / (l + x^ ir ) =T v i T /Tx for distant clusters 
(Ota & Mitsuda 2002; see section 3). 

In virialized clusters, kT vlI / pm ~ a 2 for the line-of-sight 
velocity dispersion a r of galaxies, so that 



PspecTx — T v ' lr 

is expected with 
a 2 



3 = -2^- 

A-'spcc — 9 



(6) 



(7) 



kT x ' 

where <r| al and a 2 as are the velocity dispersions of galaxies 
and gas, respectively. Equation (5) and figure 1 suggest 

I3T X =i T vir (8) 



or (3 ~ A 



spec 



for 



fvir/fc 3> 1- Otherwise, if a; V i r is 



not enough large, f3 would deviate from (3 spcc depend- 
ing on avir for given T v i T /Tx- Figure 2 shows the re- 
lation between r c /r v - 1Y (= l/x v ; r ) and j3, where we take 
T v i T /Tx = 0.6. With increasing r c /r v [ T (or decreasing 
£vir)j (3 increases monotonically. Actually, such tendency 
is found by observations (e.g., Neumann & Arnaud 1999). 

Equation (8) or (5) suggests that the Lx — Tx rela- 
tion, which is obtained from X-ray observations, would 
differ from Lx — T v ir- The L \ — T relations are compared 
in figure 3, where T — Tx in figure 3a and T — T v j r = 
PT x x 2 h ./{l + x 2 viT ) in figure 3b. 

For T — Txi we obtain the best-fit relation and 90% con- 
fidence levels L x = (l-48±g;||) x l0- 2 (kTx) 3m -° ™ with 



correlation coefficient 



xy 



0.49. On the other hand, for 



T = T V j r given by equation (5) we obtain less steep slope 

L x = (1.59±g:lf) x lO-H^vir) 2 ' 60 -" with r xy = 0.37. It 
should be noted that many of clusters deviating signif- 
icantly from the best-fit Lx — T v -„ relation are those of 
which surface brightness profiles arc far from spherical 
symmetry. The redshifts of the clusters are not taken into 
account for the analysis, because it is suggested from ob- 
servations that the redshift dependence of Lx — T is weak 
(e.g., Fairley et al. 2000). 
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Fig. 3. Lx — T relation for (a) Tx and (b) T v i r = /3Tx x^ iT / (1 + x^ iT ) . The filled circles and open squares represent the regular and 
irregular clusters, respectively (see section 3.3). Each solid line represents the best-fit relation in the form Lx = oT b with a and b 
being parameters. 
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Fig. 2. Relation between r c /r v ] r (= l/x v i r ) and /3 for 
T vil /T x = 0.6. 

3. Observational Implications 

We examine 121 clusters: 42 nearby clusters of 
z = 0.011-0.097 (Mohr, Mathiesen & Evrard 1999) and 
79 distant clusters of z = 0.10-0.82 (Ota 2001; Ota & 
Mitsuda 2002; Ota & Mitsuda 2004) with cosmological pa- 
rameters tt m = 0.3, f^A = 0.7 and Hq = 70 km s _1 Mpc -1 . 
Ota & Mitsuda (2002) mentioned the selection effect due 
to the detection limit of ASCA. Most of the clusters con- 
cerned here have rcdshifts z < 0.5 where the selection effect 
is weak; the numbers of clusters of z > 0.5 are 2 out of 50 
and 6 out of 71 for the smaller and larger core group (see 
below), respectively. 

The 121 clusters can be classified into two groups: 
"smaller core" that is smaller than 120 kpc and "larger 
core" larger than that, where 120 kpc corresponds to the 
"valley" in the distribution of core radii (figure 4). Some 



Fig. 4. Distribution of core radii of 121 clusters. One can see 
two groups of the core radii, which have peaks in 40-80 kpc 
and 140-300 kpc, respectively. 

clusters fit better to the double /3- model, i.e. superposi- 
tion of two single /3-models. For such clusters, we apply 
the double /3-model and choose the core that is dominant 
in the emissivity-averaged temperature. In other words, 
we take the component that dominates the luminosity of 
the cluster. 

In the following subsections, we discuss the possible ori- 
gin of the typical two core sizes, considering the r c — r v ; r 
relation, presence of cD galaxy, asymmetry in X-ray sur- 
face brightness and effect of radiative cooling. 

3.1. Relation between Virial Radius and Core Radius 

We assume that clusters formed by spherical collapse 
and virialized at z v - lr are observed at z Q b s = z V ir ■ We define 
r v ; r as a radius at which the average density is equal to A c 
times the critical density at z v j r , i.e. p(r v ; r ) = A c p cr it(.z v i r ) 
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Fig. 5. r c — r v j r relation. The open and filled circles repre- 
sent the clusters of smaller and larger core group, respectively. 
The thin and thick solid lines represent the best-fit relation 
(r c = ar j ) for the larger core group without statistical errors 
and that for the clusters of r c = 140—300 kpc of the Ota &i 
Mitsuda (2002) sample, respectively. The dashed line repre- 
sents r c = r v ; r /15 (see section 4). 

with A c = 187r 2 ft m (z vir ) - 427 (Nakamura & Suto 1997). 
For Tvi,. thus determined, we examine its correlation with 
r c for the 121 clusters (figure 5). 

For the larger core group, we obtain the best-fit relation 
r c = 2.95 x 10 _2 r 2 j° 5 with correlation coefficient r xy = 0.54; 
here, statistical errors are ignored, since the errors of r vir 
are not available for the Mohr, Mathiesen & Evrard (1999) 
sample. If restricted to the Ota & Mitsuda (2002) sam- 
ple for which the errors are available, we obtain r c = 

o nfi+ - 76 

(7.94™) x lO-VT; -° 54 with xVd.o.f. = 59.03/43. 

However, it should be noted that the r c distribution 
of the larger core group exhibits a tail extended to ~ 
1 Mpc (figure 4), and the r c — r v - u - relation likely de- 
pends much on the range of r c to be taken into account 
for the analysis. Actually, for the clusters in the nar- 
rower r c = 160-250 kpc and broader r c = 120-400 kpc 

q gg+0.97 

ranges, we obtain r c = (1.38±^;|f) x lO^V^ "° 56 with 

X 2 /d.o.f. = 4.03/13 and r c = (1.91±f;||) x 10~ 2 ™J-° : " 
with x 2 /d.o.f . = 26.90/32, respectively. For the intermedi- 
ate range r c = 140-300 kpc, we obtain r c = (6.31^3 8 7 24 ) x 

1 16 +1 - 21 

10" 2 r vir -° 79 with x 2 /d.o.f. = 17.25/24. From this re- 
sult, we may consider that the clusters of the larger core 
group marginally satisfy the self-similar relation (Ota & 
Mitsuda 2004). Assuming the self-similarity r c oc r vu . for 
the clusters of r c = 140-300 kpc, we obtain the best-fit 
relation r c = r V j r /16.6. 

On the other hand, for the smaller core group, we find 
no clear correlation between r c and r v ; r . This may imply 
that the origin of small cores is not explained by simple 
self-similar collapse. 



3.2. Relation between cD Galaxy and Core Radius 

cD galaxy is a giant elliptical galaxy which is located 
near the gravitational center of a cluster, and its mass 
reaches 10 11-12 M Q which is comparable to the sum of 
other member galaxies. In our sample, 50 out of 121 
clusters have central cDs or giant ellipticals. We clas- 
sify the clusters by BM-type (Bautz & Morgan 1970) 
into cD cluster (Type I and I-II) and others (Type II, 
II-III and III), using NASA/IPAC EXTRAGAL ACTIO 
DATABASE, and find 75 BM-types: 38 in nearby and 37 
in distant clusters. Although there is some question about 
the accuracy of the BM-type classification, as far as a clus- 
ter possesses one and two brightest galaxies, the cluster 
can safely be assigned to Type I and I-II, respectively (see 
e.g., Leir & Bergh 1977). 

We classify the 75 clusters further into the smaller and 
larger core groups by r c ~ 120 kpc described above. For 
each group, the fraction of clusters classified into cD clus- 
ters is shown in table la. The fraction of cD cluster is 
large in nearby clusters compared to the distant ones in 
both the smaller and larger core groups. This could be 
due to an observable limit of the surface brightness for 
the faint envelopes of distant cD galaxies. 

As for the nearby clusters, 73% of the smaller core group 
and 52% of the larger core group are classified into cD 
clusters; no distinctive difference is found between the two 
groups. The 12 (52%) Type I and I-II clusters of the larger 
core group evidently have dominant cD galaxies or giant 
ellipticals; cD galaxies are not characteristic of the smaller 
core clusters. Therefore, we may say that the origin of 
small cores are not always ascribed to the presence of cD 
galaxies. 

3.3. Relation between X-ray Symmetry and Core Radius 

The /3-model-based core size may be affected by asym- 
metry in the X-ray surface brightness. Such asymmetry 
can be due to the presence of radio lobes or cold clumps 
in the core, or caused by mergers. 

We classify the 121 clusters into "regular" and "irreg- 
ular" groups simply by the surface brightness symme- 
try. The classification of distant clusters is referred to 
Ota (2001), who examines the deviation of the surface 
brightness centroid varying the aperture radius; if the de- 
viation is greater than the 3c level, the cluster is classified 
into irregulars. For 19 nearby clusters, following Mohr et 
al. (1995) (see also Mohr, Fabricant, & Geller 1993), we 
classify clusters with wg/2 > 3a w into irregulars. For other 
20 nearby clusters, following Buote & Tsai (1996), we clas- 
sify clusters with P 1 (pk) / P (pk) > 3cr into irregulars, where 
a is estimated from the 90% confidence in their Table 2. 
Since the data are not available for three nearby clusters, 
we investigate 39 nearby and 79 distant clusters, i.e. 118 
clusters in total. 

The fractions of regular clusters are shown for the 
smaller and larger core groups in table lb. The regular 
fraction is larger in the smaller core group of the distant 
clusters, while the fraction is almost the same between the 
smaller and larger core groups for the nearby clusters. The 
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Table 1. The fractions of cD clusters and "regular" clusters. 
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Classification Nearby Distant Total 





Clusters in the smaller core group 


16/42 


(38%) 


34/79 


(43%) 


50/121 


(41%) 


(a) 


cD clusters in the smaller core group of BM . . 


11/15 


(73%) 


2/9 


(22%) 


13/24 


(54%) 




cD clusters in the larger core group of BM .... 


12/23 


(52%) 


4/28 


(14%) 


16/51 


(31%) 


(b) 


Regular clusters in the smaller core group 


10/15 


(67%) 


28/34 


(82%) 


38/49 


(78%) 




Regular clusters in the larger core group 


15/24 


(63%) 


17/45 


(38%) 


32/69 


(46%) 



(a) 75 clusters are classified into the smaller and larger core groups. 

(b) 118 clusters are classified into the smaller and larger core groups. 



difference between the distant and nearby clusters may be 
due to mergers or be interpreted as that distant clusters 
are well relaxed compared to nearby ones, though there 
could be an observational effect such that detection of sub- 
structures in distant clusters are difficult. As for the core 
sizes, however, no clear difference is found in the nearby 
clusters. Therefore, it is unlikely that X-ray asymmetry 
is characteristic to small core clusters. 

3.4. Relation between Radiative Cooling and Core 
Radius 

In the cores of some clusters, the gas undergoes signifi- 
cant radiative cooling, and then inflow of the ambient gas 
may occur to compensate the pressure loss inside (e.g., 
Fabian 1994). This process forms a peaked density profile 
responsible for the smaller cores. 

The cooling time at the cluster center due to 
bremsstrahlung is written as 



^COol ; (9) 

qsn e o 

where n e o is the electron number density and qg is a co- 
efficient for the emissivity q^^/T. In figure 6a is shown 
icooi against the core radius r c for the 121 clusters. The 
best-fit relation is obtained as £ coo i = 2.24 x 10 3 r 325 with 



with 



= 0.44 for the smaller core group and t coo i = 5.50 r^ 07 



xy 



0.55 for the larger core group. 



In Section 2, we argue that Tx is not always a good 
quantity to represent T vlr in the context of the /3-modcl. 
So that, we define another cooling time 



t 



cool 



qsn e o 



1 



1/2 



t, 



COOl 3 



(10) 



which depends on [3. The cooling time t' cool is shown in 
figure 6b. We obtain the best-fit relation t' cooi = 3.80 x 
10 2 r 2 71 with r xy = 0.51 for the smaller core and t' cool = 
5.01 r\ w with r xy = 0.75 for the larger core groups. The 
relation of f/ ool with r c is clearer than that of t coo i for the 
larger core clusters which have long cooling time. 

In comparison with the larger core group, the relation of 
r c with cooling time is not clear in the smaller core group. 
As discussed in section 3.1, neither the self-similar relation 
between r c and r v j r is clear in the smaller core group. 



These may suggest that the clusters of small cores are 
affected by radiative cooling after the self-similar collapse. 

Also cDs or galaxies in the central region can alter the 
core structure through their AGN activities, supernovae 
feedback or galactic winds. However, as discussed in sec- 
tion 3.2, the difference in the relation with such galaxies 
is not clear between the smaller and larger core groups as 
compared to the relation with the cooling time. Moreover, 
heat or pressure due to such activities of central galaxies 
likely works against the core getting smaller, while the 
relevant metal enrichment enhances radiative cooling. 

If the gas-mass fraction in the cooling core is nearly 
constant among the clusters, a scaling law follows the self- 
similarity r v ; r cx r c , as (see section 2 and equation (18)) 



1/2 -1 
vir n e0 



OCr e r vir 0Cr c- 



(11) 



This dependence is seen in figure 6b for the clusters of 
the larger core group. For large core clusters of t' cool > 
in, radiative cooling is not yet significant and their cores 
remain to follow the self-similar relation. Otherwise, if 
cooling is significant, the core would evolve to deviate from 
the self-similarity This is likely the case for the smaller 
core group: no clear correlations between r c and ?* v i r in 
figure 5. 

In order to see the effect of cooling on the smaller core 
group clusters, we define the cooling radius based on the 
[3- mo del, as 



^cool ^cool,10 

or 



-2/3/3 



cool 



cool, 10 



-2/3/3 



-1 



- 1 



1/2 



1/2 



(12) 



(13) 



for t 



cool. 10 



or t' 



cool, 10 



in units of 10 Gyr, and investigate 
70 clusters of t CO oi,io < 1- The cooling radius is plotted 
against the core radius in figure 7. No clear correlations 
are seen in r coo i - r c or r' cool - r c . 

If t' cool is significantly shorter than in, we have an ap- 
proximate relation 

'cool 01 'c'cool 01 'c 

where, in the last expression, t' 



(14) 

cx r 2 - 7 suggested from 
figure 6b is taken into account, though the correlation is 
quite weak. Consequently, for various values of j3 from 
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Fig. 6. Relation of r c with (a) t coo i and with (b) t' ^. The solid lines represent the best-fit relation, t coo l (or t' j) = ar', for 
the smaller core group (open circles) and for the larger core group (filled circles). Statistical errors are not taken into account for 
the analysis, since the errors of cooling time are not available for the Mohr, Mathiesen & Evrard (1999) sample. The dashed line 
represents the Hubble time tjj- 



Table 2. Resolutions h and systematic errors 0(h 2 ). 



Asph 


ha.oi* 


/l0.05* 


ho.i* 


0(h") 




(kpc) 


(kpc) 


(kpc) 


(%) 


30000 


31.5 


18.4 


14.6 


0.961 


50000 


26.6 


15.6 


12.3 


0.705 



* The gas density n g (cm 3 ). 

~ 0.6 to ~ 1 (figure 1), r c dependence of r' cool would be 
smeared and the self-similar relation might be lost for the 
smaller core group. This suggests that cooling in the core 
is likely connected with the origin of small cores of clusters 
which are out of self-similarity. 

4. Hydrodynamical Approach 

Using an SPH (Smoothed Particle Hydrodynamics) 
code based on Hcrnquist & Kats (1989) and 
Monaghan (1992), we simulate the /3-model and in- 
vestigate the behavior of the isothermal intracluster gas 
in order to make sure the arguments in the previous 
sections. Our main concern here is the self-similar 
relation r c oc r v j r and the effect of the potential due to 
cD galaxy on the isothermal gas without cooling; the 
effect of radiative cooling will be discussed in a separate 
paper. Table 2 shows the resolution and systematic error 
estimated for the calculations. 

We consider that a cluster consists of galaxies, gas and 
dark matter with the mass ratio M ga i : M gas : A/dm = 1 : 
5 : 30. The Momentum equation for the gas is written as 



d,Vj 1 

-1 = VPi + aJ 

at p t 



V$i - Gi 



(15) 



where the first and third terms on the right hand side 
represent the force due to the pressure gradient and 



self-gravity, respectively. The second term is the artifi- 
cial viscosity, for which we adopt the Monaghan- Gingold 
value (Monaghan & Gingold 1983; see also Hernquist & 
Kats 1989). The last term represents the force due to the 
components other than the gas, i.e. the dark matter and 
member galaxies, 



-Gi 



GM cx (n) n 
r? I n 



(16) 



where M ex (rj) is the mass of the components within radius 
r%. For their distributions, we take an approximate King 
profile 

-3/2 



Pk(t) = pko[1 + (r/r K ) 



(17) 



with pYi(r > r v j r ) = 0, where pxo is given for each compo- 
nent according to the mass ratio and tk is taken to be the 
initial core size of gas. 

For the initial distribution of the gas we set the /3-model 
with /3=1, rigo = 0.04 cm~ 3 and r v i T /r c = x v i T = 15 (see 
figure 5), and calculate the evolution to be hydrostatic 
with the gas being isothermal. The initial value of the gas 
temperature T is taken to be T V - 1T , considering (3x^/(1 + 
1 at the initial state. 



< ir ) = 0.996 



4-1- Comparison with the (3-model 

We carry out calculations for four clusters of the virial 
radius 2.3, 3.0, 3.6 and 4.0 Mpc until 3.0 Gyr, which is 
enough longer than the free fall time ~ 1.2 Gyr for n g o = 
0.04 cm" 3 , and examine the resultant intracluster gas. 
We find that the gas profile is well reproduced by the (3- 
model with the parameter values shown in table 3. The 
resultant (3 is smaller than unity (initial value), yet the 
relation T v j r ~ (3T is attained for the gas temperature T, 
as expected from the argument in section 2. 

Also the self-similar relation between r c and r v j r is kept 
in spite of that (3 becomes smaller than unity, as shown in 
figure 8a. This implies that the core radius well reflects 
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Fig. 7. Relation of r c with (a) r coo i and with (b) r' >. For a double /3-model cluster, each core radius is plotted: the filled circle 
and open square represent the dominant core and another (minor) one, respectively. The single /3-model clusters are represented by 
the filled circles. 



Table 3. Best-fit parameters of the /3-model. 





Model 




After 3.0Gyr 


M vir 


V yir 


V 




(1O 15 M ) 


(Mpc) 


(kpc) 


(kpc) (cm -3 ) 


0.48 


2.30 


150 


150 2.00 0.83 


1.03 


3.00 


200 


200 2.04 0.86 


1.72 


3.60 


240 


240 2.03 0.87 


2.50 


4.00 


270 


270 2.10 0.88 



the core of the gravitational potential consistently with 
the /3-model. The consistency can be examined in another 
way as follows. When the gas follows the /3-model, the gas 
density at the cluster center is written as 

9fc/3T/ s o 



(18) 



where f g o is the gas-mass fraction at the cluster center. 
If the self-similar relation r c attained, T V i r /r^ ~ 

fiT/r^. becomes independent of r c from equation (4) and 

3 



M vir oc r. 3 



lr oc , which is obtained also from the calcula- 



tions (see table 3). Accordingly, n g0 /f g0 should be inde- 
pendent of r c , if the calculated structure is consistent with 
the /3-model. This behavior is seen in figure 8b, where 
n g o/f g o is plotted against r c for the four calculations. 

4-2. cD galaxy 

Intending to estimate the maximal effect of a central cD 
galaxy, we take a large cD galaxy with the core radius ?'k = 
20 kpc into account for the cluster potential. We consider 
a cluster of r v - K = 3 Mpc, r c = 200 kpc and n g o = 0.04 cm -3 
and simply replace the galaxy component by the cD one 
so as for r v j r , M vir , T vir and the gas-mass fraction to be 
the same values as those in the case without the cD for the 
same mass ratio; this also allows to keep the resolution for 
the computation resources. Although the potential due to 



other member galaxies is ignored here, we make sure that 
their contribution is negligibly small compared to other 
components. The integrated mass of each component is 
shown in figure 9. 

Figure 10 shows the gas profile after 1.2 Gyr with N = 
50000 SPH-particles. The density increases toward the 
cluster center in comparison with the initial profile. The 
best-fit /3-model (see table 4) to the resultant profile gives 
r c = 150 kpc smaller than the initial core size, while the 
profile at > 100 kpc varies little from the initial one. This 
means that the presence of a central cD galaxy appears 
merely as an excess within 100 kpc and little affects the 
outer region. Therefore, we try to fit the gas profile using 
the double /3-model as follows. 

We first obtain the best-fit parameters n g o, r c and (3\ 
for the outer component, r = 0.1-2.0 Mpc, and then in- 
clude the inner component to obtain ri g o,2 and r c ^ with 
n g o, r c and /3i fixed, where we assume 02 = Pi as done 
in the double /3-model analyses by Mohr, Mathiesen & 
Evrard (1999) and Ota & Mitsuda (2002). The inner core 
radius thus obtained, r c , 2 = 41 kpc, is about twice the core 
radius tk = 20 kpc of the cD galaxy. The extended core 
of the gas compared to the cD potential is likely due to 
the pressure gradient. As seen in figure 4, the core radius 
~ 40 kpc is comparable to the smaller end of the smaller 
core group, but is too small to explain the distribution 
40-80 kpc because the typical core size of cD galaxies is 
less than ~ 20 kpc (see e.g., Sparke & Gallagher 2000). 

5. Concluding Remarks 

We summarize the points discussed in the present work. 

1. The isothermal X-ray temperature of a cluster, 
which is inferred from the X-ray surface brightness 
in the context of the /3-model, does not always rep- 
resent the virial temperature of the intracluster gas, 
even if the cluster is well relaxed. The virial tem- 
perature T vir is given by /3T x ^ ir /(l + ^.J ~ P T x 
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Fig. 8. (a) r c — r v ; r and (b) n g g/f g Q — r c relation. The dotted lines represent the initial values: (a) r c = r v ; r /15 and (b) 
n g0 /fg = 0.04/0.139 cm -3 . The circles represent the results after 3.0 Gyr with N = 30000 SPH-particles. 



Table 4. Best-fit parameters of the single and double /3-models. 





Sinj 


de /?-model 


Double /?-model a 




V 


n g0 Pi 


r c ,2 


%0,2 @2 


Model 


(kpc) 


(cm -3 ) 


(kpc) 


(cm -3 ) 


I/? 


150 


5.27 0.91 






1/3 (MOO kpc) ... 


181 


4.04 0.95 






2/3 






41 


2.66 



a Parameters of the outer component are fixed at the values obtained 
outside 100 kpc. 




Fig. 9. Initial profile of the integrated mass: dark matter 
(dashed line), gas (dot-dashed line), cD galaxy (dotted line), 
and total (solid line). 



Fig. 10. Gas density profile: the circles represent the profile 
after 1.2 Gyr with TV = 50000 SPH-particlcs and the bars 
attached represent the errors (smoothing lengths). The lines 
represent the best-fit single and double /3-models. 
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(for Xvir 3* 1) rather than Tx itself, where Tx and (3 
are the isothermal gas temperature and profile-slope 
parameter, respectively, and x v i T is the ratio of the 
virial to core radii, i.e. x vlI = r v - n -/r c . 

2. For 121 clusters observed by ROSAT and ASCA, the 
luminosity-temperature relation Lx — f3Tx , which 
should represent Lx — Tvir, is less steep than Lx ~ 
Tx- This result may imply that the /3-correction 
to Tx is important in the study of the properties 
expected from the self-similar collapse. 

3. The 121 clusters can be classified by the core 
size into two groups, as suggested by Ota & 
Mitsuda (2002). No clear correlation is seen be- 
tween r c and r v ; r for the small core (~ 40-80 kpc) 
clusters, while the self-similarity r c oc r v ; r is attained 
marginally for the larger core (~ 140-300 kpc) clus- 
ters. Neither clear relation with cD galaxy is found 
for the small core clusters. 

4. It is found, however, that the small core clusters 
have cooling time significantly shorter than the 
Hubble time, while the larger core clusters have 
longer cooling time. These suggest that radiative 
cooling in the core may be responsible for the small 
cores or breaking the self-similarity r c oc r v ; r , which 
is expected from the collapse and seen in the larger 
core clusters. 

5. Hydrodynamical calculations of the isothermal gas 
in clusters, which consist of the dark matter, galax- 
ies and gas, show the self-similar relation r c oc r V i r 
and T v j r ~ (3T consistently with our analysis of the 
(3- model (section 2). Calculations also show that a 
large cD galaxy located in the cluster center affects 
the gas profile within the region ~ 100 kpc. Thus, 
under the presence of a central cD galaxy, the intra- 
cluster gas may be represented better by the double 
/3-model than the single one. 

6. However, the inner core (~ 40 kpc) due to the cD 
galaxy is too small to account for the core-size dis- 
tribution of the smaller core group. Finally, from 
the analysis of the 121 clusters as well, it is unlikely 
that the origin of the smaller cores is explained by 
the presence of a central cD galaxy. 

The authors would like to thank Naomi Ota for her 
making the cluster data available and helpful suggestions 
for their analysis. 
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